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Abstract 

' Context. Low and intermediate mass stars lose a significant fraction of their mass through a dust-driven wind during the Asymptotic Giant 
, Branch (AGB) phase. Recent studies show that winds from late-type stars are far from being smooth. Mass-loss variations occur on different 
• time scales, from years to tens of thousands of years. The variations appear to be particularly prominent towards the end of the AGB evolution. 
The occurrence, amplitude and time scale of these variations are still not well understood. 

Aims. The goal of our study is to gain insight into the structure of the circumstellar envelope (CSE) of WX Psc and map the possible variability 
of the late-AGB mass-loss phenomenon. 

Methods. We have performed an in-depth analysis of the extreme infrared AGB star WX Psc by modeling (1) the CO J = 1-0 through 7-6 
rotational line profiles and the full spectral energy distribution (SED) ranging from 0.7 to 1300 yum. We hence are able to trace a geometrically 
extended region of the CSE. 

Results. Both mass-loss diagnostics bear evidence of the occurrence of mass-loss modulations during the last ~2000 yr In particular, WX Psc 
went through a high mass-loss phase (M~ 5 10"^ Mo/yr) some 800 yr ago. This phase lasted about 600 yr and was followed by a long period of 
low mass loss (M~ 5 10"* Mg/yr). The present day mass-loss rate is estimated to be ~ 6 10"* Mg/yr 

Conclusions. The AGB star WX Psc has undergone strong mass-loss rate variability on a time scale of several hundred years during the last 
few thousand years. These variations are traced in the strength and profile of the CO rotational lines and in the SED. We have consistently 
I simulated the behaviour of both tracers using radiative transfer codes that allow for non-constant mass-loss rates. 

Key words. Line: profiles. Radiative transfer. Stars: AGB and post-AGB, (Stars): circumstellar matter. Stars: mass loss. Stars: individual: WX 
Psc 



1. Introduction 

The interstellar medium (ISM) is slowly enriched in heavy el- 
ements through material ejected by evolved stars. Such stars 
lose mass through a stellar wind, which is slow and dusty for 
cool giants and supergiants, or through supernova explosions. 
The ejected material merges with the interstellar matter and is 
later incorporated into new generations of stars and planets (the 
so-called 'cosmic cycle'). The gas and solid dust particles pro- 
duced by evolved stars play a major role in the chemistry and 
energy balance of the ISM, and in the process of star formation. 

Low and intermediate mass stars (M ~ 0.6 - 8 Mq) are par- 
ticularly important in this process since their mass loss dom- 
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inates the total ISM mass. About half of the heavy elements, 
recycled by st ars, originate from stars of this mass interval 
( Maedeilll992 ). Moreover, these stars are the dominant source 
of dust in the ISM. When these stars leave the main sequence, 
they ascend the red giant and asymptotic giant branches (RGB 
and AGB). During these evolutionary phases, very lar ge ampli- 
tude and long period pulsations (e.g. i BowenI 1988 ) lift pho- 
tospheric material to great heights. In these cool and dense 
molecular layers dust grains condense. Radiation pressure ef- 
ficiently accelerates these dust grains outwards and the gas is 
dragged along. The resulting mass-loss rate is much larger than 
the rate at which the material is burned in the core and domi- 
nates the stellar evolution during this phase. Despite the impor- 
tance of the mass-loss process in terminating the AGB stellar 
evolutionary phase and in replenishing the ISM with newly- 
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produced elements, the nature of the wind dr iving mechanism 
is still not well understood (e.g. IWoitkell2006h . 



Recent CO and scattered light observations of AGB objects 
sho w that winds from late-tvpe stars are far from being smooth 
(e.g. Mauron & Huggins 200a). Density variations in the cir- 
cumstellar envelope (CSE) occur on different time scales, rang- 
ing from years to tens of thousands of years. Stellar pulsations 
may cause density oscillations on a time scale of a few hundred 
days. A nuclear thermal pulse may be responsible for variations 
every ten thousand to hundred thousand years. Oscillations on 
a time scale of a few hundred years are possibly linked to the 
hydrodynamical properties of the CSE (see also Sect. 15.3b . 

The mass-loss variations appear to be particularly more 
prominent when these stars ascend the AGB and become more 
luminous. It is often postulated and in rare cases observed (e.g. 
Justtanont et~aDl 19961: Ivan Loon et aDl2003i) that the AGB evo- 
lution ends i n a very high mass-los s phase, the so-called super- 
wind phase dlben & Renzinilll983b . The superwind mass-loss 
determines the quick ejection of the whole envelope and the 
termination of the AGB phase. 

To enlarge our insight in this superwind mass-loss phase, 
and in possible mass-loss variations occurring during this 
phase, we focus this research paper on the study of WX Psc 
(=IRC-Hl0011, IRAS 01037-H1219, CIT 3), one of the mos t 
extreme infrared (IR) AGB objects known ( Ulrich et alJll966l) . 
WX Psc belongs to t he very late-type AGB stars having a spec- 
tral type of M9- 10 ( Lockwoodll985 ). Given this spectral type 
and its very red color (V-K 18), its effective tem perature can 
be estimated to be <, 2500 K (iHofmann et al.ll200lh . WX Psc is 
an ox ygen-rich, long-period variable (P ~660days; iLe Bertre 
1993h . 

The Infrared Space Observatory (ISO) data show that WX 
Psc is surrounded by an optically thick dust y shell, of which 
the presence of crystalline silicate features ('Suh"2002') is in- 
dicative of a high mass-loss rate (Kemceret al. 2001). In the 
circumstellar envelope ( CSE) SiO (iDesmurs et alJl2000l) . H2O 
dBowers & HagenI 1 1984 ) and OH (" Olnon et all 1198 0*) maser 
emission lines are detected with expansion velocities ranging 
from 18 to 23kms"'. 

Estimates of the mass-loss rate, based on the modeling of 
CO rotational line intensities or (near)-infrared (dust) observa- 
tions, span a wide range. Derived values lie between 1.4 x 10"^ 
and 1.3 X 10""* Mq yr"' (see Table[T]). The quoted gas mass-loss 
rates are however often based on one diagnostics (e.g., the in- 
tegrated intensity of one rotational transition of CO) and only 
trace a spatially limited region. To clarify our understanding 
on the mass-loss process, it is crucial to study different trac- 
ers of mass loss covering a geometrically extended part of the 
CSE. We therefore have performed an in-depth analysis of the 
CO J — 1—0 through 7—6 rotational line profiles in combination 
with a modeling of the spectral energy distribution (SED). We 
demonstrate that this combined analysis gives us the possibility 
to pinpoint the (variable) mass loss of WX Psc. 

In Sect. |2] we describe the data sets used in this paper. The 
method of analysis is outlined in Sect. [3] and the results from 
both the CO line profile modeling and SED modeling are pre- 
sented in Sect. |4] The results are discussed in Sect. |5] and are 



compared with results obtained in other studies. The conclu- 
sions are formulated in Sect. |6l 

2. Observations 

The analysis is based on ( i) rotational CO line profiles found in 
the literature and ( ii) the SED in the range from 0.7 to 1 300 fim. 
The applied CO-dataset is indicated with the superscript t in 
Table [T] Main emph asis is put on the high-quality data ob- 
tained bv lKemper et al. (2003) on the 15 m JCMT telescope on 
Mauna Kea, Hawai. The half power beam width of the different 
receivers was 19.7" for CO (2-1), 13.2" for CO (3-2), 10.8" 
for CO (4-3), 8.0" for CO (6-5), and 6.0" for CO (7-6). The 
absolute flux calibration accuracy was estimated to be 10 % for 
the CO (2-1), (3-2) and (4-3), and -30 % for the CO (6-5) 
and (7-6) line because the flux standards, that are available for 
these higher excitation lines, are less reliable. We note that in 
principal pointing uncertainties, beam effects, atmospheric see- 
ing, and systematics at these frequencies can conspire against 
you to produce an error that is larger than 30 %. Unfortunately 
we don't have enough statistics information on standards etc. 
to know for sure how bad a problem this was for the RxE and 
RxW/D receivers at the time of the observations. 

The SED of WX Psc is that of a strongly reddened Mira- 
type variable star with a period of 66 days and a K band am- 
plitude of 3.4 mag ( Samus et al ] E004l) . Two full spectral scans 
from both spectrometers (Short Wavelength Spectrometer 
{SWS) and Long Wavelength Spectrometer (LWS)) onboard the 
Infrared Space Observatory (ISO) satellite are available. Due 
to the pulsation of the star the flux levels of the SWS spec- 
tra differ by a factor ~2. We opt to use the SWS spectrum 
taken on 15 Dec 1996 (tdt 39502217) and the LWS spec- 
trum from 15 Jun 1997 (tdt 57701103). Although these spec- 
tra have not been obtained simultaneously, the flux levels in 
the overlapping region between SWS and LWS agree reason- 
ably well because the observing dates correspond to phases 
0.38 and 0.66, i.e. roughly symmetric compared to the mini- 
mum at 0.5. We have further used th e 0.97-2.5 jt/m near-IR 
spectrum from lLancon & Wood ( 2000l) scaled to the flux-level 
of the SWS spectrum. The IR spectroscopic data have been 
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suppl e mented by optical and IR photometric dat a jCutri et al 
I2OO3I; iDvck et all fl974 
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sub- 

mm photometry (Sookaetal. 1989) and radio measurements 
( Walmsley et al.i.l991:.Dehaes et al..2006i) in order to construct 
the fuU SED (see Fig. [B. 

In the wavelength region up to ~8 pm, photospheric ab- 
sorption features and emission features arising from density 
enhancements in the outer atmosphere caused by pulsational 
levitation are visible. Specifically, the CO first overtone and 
H2O vi symmetric stretch and V3 asymmetric stretch band fea- 
tures are arising between ~ 2 and 3.8//m, the small absorp- 
tion peak around 4.05 pva is caused by the SiO first overtone, 
with the stronger absorption peak at slightly larger wavelengths 
(at ~5 pm) mainly being due to the CO fundamental vibration- 
rotation band. From 6.21 pva onwards, the H2O V2 vibrational 
bending band is the main absorption feature. However, for the 
purpose of determining the mass-loss history of WX Psc, a de- 
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Figure 1. SED of WX Psc (black), together 
with the best fit dust model (gray) and the 
best fitting constant mass-loss rate model 
(dashed). The inset shows the same SEDs 
using logarithmic axes in order to empha- 
size the far-IR slope. Note that the large dis- 
crepancies between the observed and mod- 
eled SEDs between l-8yum are due to the 
representation of the stellar atmosphere by 
a blackbody. The wavelength range which 
is dominated by (unmodeled) photospheric 
features is to the left of the dotted line at 
8yum. The main absorption and emission 
features in the 1-8 yum range are due to 
H2O, SiO and CO. The dust absorption and 
emission is best traced beyond 8//m (see 
text for more details). The dotted line in the 
inset gives an indication of the level of the 
photospheric contribution to the total flux. 



Table 1. Mass-loss rates for WX Psc as found in literature. Data used in this paper are indicated with the superscript 
appropriate telescope diameter (in meter) and beam size (in ") are listed in the last column. 
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Observ. data 


M [Moyr"'] 


Comments 


Ref. 


Diam. - beam 


CO (1-0) 


1.2 X 10^^ 


from analytical expression 


Knapp & Morris (1985) 7m - 


100" 


CO (1-0) 


2.4 X 10-5 


theor. model of Morris (1980) 


Sopkaetal. (1989) 


20m 


-42" 


CO (1-0) 






Mareulis et al. (1990) 


14m 


-45" 


CO (1-0) 






Nymanetal. (1992)+ 


20m 


-42" 


CO (1-0) 


8.5(±2.39) X 10-'' 


from analytical expression 


Loup etal. (1993) 






CO (1-0) 


1.0 X 10-^ 


self-consistent T(r) + NLTE radiat. transfer Justtanont et al. (1994) 






CO (2-1) 


2.6 X 10-^ 


from analytical expression 


Knapp et al. (1982) 


10m 


-30" 


CO (2-1) 


1.0 X 10-5 


large dilference between data and model 


Justtanont et al. (1994) 






CO (2-1) 


5.0 X 10-5 


from analytical expression 


Jura (1983) 






CO (2-l),(3-2), 


0.14- 8.0 X 10-« 


from model Justtanont et al. (1994) 


Kemper et ak (2003)+ 


15m 


- (see text) 


(4-3),(6-5),(7-6) 












L- [12yum] 


2.0 X 10-5 


dust modeling 


Kemper et al. (2002) 






OH maser 


2 X 10-5 


usins Eq. (3) in Bowers et al. (1983) 


Bowers etal. (1983) 






J, H, K' + 


1.3 - 2.1 X 10-5 


spherical dusty wind with 


Hofmannetal. (2001) 






SED + ISI(ll/im) 




p cx f--' 5 from 137 Ri, on 








J, H, K' + 


1.3 X lO-'' 


spherical dusty wind 


Tevousian et al. (2004) 
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y, H, K' + 


9x 10-*^ last 550 yr 


spherical dusty wind -1- cocoon 


Vinkovic et al. (2004) 
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2 X 10-5, _ 5500 


yr 









tailed analysis of the atmospheric behaviour of this target is 
not a prime issue. In the framework of the dust modeling pre- 
sented here, we simply represent the central star by a black- 
body. The assumption of a blackbody for the underlying star 
does not influence the conclusions since dust emission and ab- 
sorption dominates over photospheric effects beyond 8 jum (see 
Fig.[T]i. We focus on several dust characteristics of the SED and 
IR spectroscopy to estimate the mass loss. The most significant 
of these are 7 ) the average reddening of the star; 2 ) the self- 



absorption of the lOjum silicate feature; 3) the relative strength 
of the 10 and 20jum bands, and 4) the slope in the mid- to far- 
IR. These observables trace the relative amounts of dust near 
and further away from the star, i.e. the columns of warm and 
cooler material along the line of sight. 

3. Analysis 

To constrain the dust-driven wind structure from its base out 
to several thousands of stellar radii, we use both molecular 
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line fitting of the CO rotational J - I - through 7-6 line 
profiles and dust radiative transfer simulations of the global 
SED and IR spectrum. The spectral lines are modeled us- 
ing GASTRoNOoM, the dust continuum using modust. We 
briefly discuss both programs in Sect. 13.11 Special emphasis 
is given on the model selection procedure and on the measure 
of the goodness-of-fit in Sect. 13.1.31 The results are presented 
in Sect.m 

3.1. Theoretical code 
3.1.1. GASTRoNOoM 

The observed line profiles provide information on the thermo- 
dynamical structure of the outflow of WX Psc. For a proper 
interpretation of the full line profiles, we have developed a the- 
oretical model (GASTRoNOoM - Gas Theoretical Research 
on Non-LTE Opacities of Molecules) which first (1) calcu- 
lates the kinetic temperature in the shell by solving the equa- 
tions of motion of gas and dust and the energy balance si- 
multaneously, then (2) solves the radiative transfer e quation 
the co-moving frame (CMF, iMihalas et alJ Il975h using 



m 



co-moving frame (CMF, 
the Approximate New ton-Raphson operator as developed by 
Schonberg & Hempe and computes the NLTE level- 

populations consistently, and finally (3) determines the observ- 
able line profil e by ray- trac ing. A full description of the code 
is presented in iDecin et al.l (i2006.) . The main assumption is a 
spherically symmetric mass loss. The mass-loss rate is allowed 
to vary with radial distance from the star. Note that the gas and 
drift velocity are not changed accordingly. Adopting a value 
for the gas mass-loss rate, the dust-to-gas ratio is adjusted until 
the observed terminal velocity is obtained. For details on the 
thermal balance, the gas and dust velocity, and the t reatment 
of variable mass loss we r efer to lOecin et al.l (l2006t) . We al- 
ready note that for solving the radiative transfer equation in the 
CMF-frame in the case of WX Psc, we have used a value for 
the mass-loss rate, M, equal to 1 x 10"^ M^lyx. 



3.1.2. MODUST 

We have used the code modust (' Bouwman et al.l2000 ^ to model 
the transfer of radiation through the dusty outflow of WX Psc. 
This code has been applied extensively to model not only 
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Dijkstra et al] 120031: 



Diikstra et al. 2006), but als o solid state par- 
ticles in proto-planetary disks and comets (e.g. Bou wman et al 



20031) . In short, the code assumes a spherical distribution, such 
as an outflow or shell, of dust grains in radiative equilibrium 
that are irradiated by a central star. Here, we represent the pho- 
tospheric energy distribution of WX Psc by a blackbody (see 
Sect.|2]). A special circumstance in the modeling of WX Psc is 
that we account for two distinct dust shells (see Sect. 14.2b . For 
the inner shell we do a consistent radiation transfer solution, 
for the outer shell we assume that the dust medium is optically 
thin and that it is irradiated by the light emerging from the in- 
ner shell. We have verified that indeed the outer shell of WX 
Psc is optically thin. 



3.1 .3. Parameter estimation, model selection, and the 
assessment of goodness-of-fit 

The observational line data are subject to two kinds of uncer- 
tainties: (1) the systematic or absolute eiTors, with variance 
cr~^^, arising e.g. from uncertainties in the correction for vari- 
ations in the atmospheric conditions, and (2) the statistical or 
random errors e, with variance cr^taf which reflect the variabil- 
ity of the points within a certain bin. Both kind of errors have 
to be taken into account when estimating the model parame- 
ters. In case of the CO line profiles, the statistical errors (rms) 
are smaller than the systematic errors. We hence will construct 
a statistical method giving a greater weight to the line profile 
than to the integrated line intensity. Moreover, as demonstrated 
in Decin et al. (2006), especially the line profiles yield strong 
diagnostics for the determination of the mass-loss history. 

To determine the mass-loss history of WX Psc from the 
CO line profiles, a grid of ~300000 models was constructed. 
In a first step, the models for which the predicted line profiles 
did not fulfill the absolute-flux calibration uncertainty criterion 
were excluded. This absolute flux criterion does not only in- 
clude the absolute flux uncertainty on the data as specified in 
Sect. |2] but also incorporates the noise on the data. From then 
on, the other models were treated equally when judging the 
quality of the line profile shape. To do so, a scaling factor based 
on the ratio of the integrated intensity of the observed data to 
the integrated intensity of the predicted data was introduced to 
re-scale the observational data. From a statistical standpoint, 
this scaling factor merely is an additional model parameter, ac- 
counting for the calibration uncertainties. 

Assume a model for the re-scaled observed spectrum y, at 
velocity i, i - 1,2, . . . ,n with expected mean value E(3',) = fii, 
representing the theoretical spectrum of the target. The statisti- 
cal measurement errors are assumed to be normally distributed 
with mean and variance taken to be constant at all frequency 
points i, i.e. (Tstat; = CTstat- Hence, the model is assumed to fol- 
low a normal distribution 



yi = /Ji + Ei, where ~ A^(0, a^^J. 



(1) 



A key tool in model selection and fitting, and in the assessment 
of goodness-of-fit, is the log-likelihood function which, for the 
normally distributed case, takes the form {, defined as: 



-ln(o-st,t,)-ln(V2^)- - 



1 (yp-f^jJ 



(2) 



j=\ \i=\ 

where Mines represents the total number of rotational lines be- 
ing 6, a nd n , the number of frequency points in the y-th line 
profile dCox & HinklevI Il990h . The model which best fits the 
observed data maximizes the log-likelihood distribution. 

One can derive three p roperties from the log-likehood dis- 
tribution (see Welshl[l996l) . First, for a given parametric model, 
maximizing t leads to the so-called maximum likelihood esti- 
mator for the model parameters. Denote by €,„ the value of 
the log-likelihood function at maximum. Standard eiTors are 
deduced as the square root of the diagonal of the negative in- 
verse information matrix. The information matrix is the matrix 
of second derivatives of the log-likelihood function w.r.t. the 
model parameters. 
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Table 2. (1 - a) 100% quantile of the ;f p-distribution, with p 
degrees of freedom and a significance level a - 0.05. 



p 




P 




1 


3.8414588 


11 


19.675138 


2 


5.9914645 


12 


21.026070 


3 


7.8147279 


13 


22.362032 


4 


9.4877290 


14 


23.684791 


5 


11.070498 


15 


24.995790 


6 


12.591587 


16 


26.296228 


7 


14.067140 


17 


27.587112 


8 


15.507313 


18 


28.869299 


9 


16.918978 


19 


30.143527 


10 


18.307038 


20 


31.410433 




Second, the log-likelihood function can be used to deter- 
mine (1 - a) 100 % confidence intervals for the model param- 
eters, where a is the significance level. Precisely this is done 
through the method of profile likelihood, and the correspond- 
ing intervals are termed 'profile likelihood confidence inter- 
vals'. Let us focus on one of the parameters, (f say. The pro- 
file hkelihood interval is defined as the values (f for which 

= ({ip,6{(f)) is sufficiently close to i,„. We use 6{(f) as the 
parameter that maximizes the log-likelihood function subject 
to the constraint that if is kept fixed. Sufficiently close is to be 
interpreted in the precise sense that 2(^„, - i^p) < x^(<^), where 
xlicf) is the (1 - a) 100% quantile of the;t'j distribution. The 
same procedure can be used for a group of p parameters si- 
multaneously. In this case ip is /^-dimensional and the quantile 
is now of the ;t'^(a')-distribution (see Table|2|. In our case, the 
confidence intervals are influenced by the grid spacing, but they 
converge to the true intervals if the grid spacing becomes arbi- 
trarily fine. 

Third, the log-likelihood function can be used to compare 
two different models, with a different number of parameters, p 
and p + r, as long as they are nested. The more elaborate model 
is preferred if 2{tf„^'^ - i'^,) > x^icf), otherwise preference points 
towards the simpler model. 

One should realize that the outlined log-likelihood method 
does not force the accordance between observations and pre- 
dictions to be better than a certain pre-specified level in or- 
der to be acceptable. This is in contrast with the often used 
reduced-;^'^-method, which requires that for a good fit the value 
of the Teduced-x^ should be lower than 1 (or 2). Using the log- 
likelihood method, one is selecting the theoretical predictions 
(and corresponding confidence intervals for the model parame- 
ters) with the best fit to the data in the framework of the applied 
theoretical model. 

As an example, we describe the selection of the best model 
and the determination of the 95 % confidence intervals of the 
stellar and mass-loss parameters as derived from the CO line 
fitting procedure (see Sect. 14. lb . In Fig.|2]an example is shown 
of a mass-loss profile with 4 M-phases and 3 small interfaces 
(each described by 2 parameters (/?,-, Mi)). This mass-loss pro- 
file is hence described by 6 pairs (/?,, M,), indicated with dia- 
monds in Fig. 121 Other free parameters in the grid computation 



Figure 2. Upper: Estimated temperature profile, middle: esti- 
mated velocity structure, and bottom: estimated mass-loss rate 
M{r) for WX Psc as a function of radial distance to the star 
(black line). The shading region in the bottom panel displays 
the confidence intervals of the estimated mass-loss parameters. 

were the stellar temperature T*, dust condensation radius Thinner, 
and outer radius of the CSE 7?outer, resulting in 15 free param- 
eters in total. We note that for the modelling of WX Psc the 
value of M at Thinner was specified to stay constant till the first 
decrease in mass loss at ~3QR-^, and also from 15507?* on the 
mass-loss was kept constant until /?outer- This is however not a 
prerequisite to the code, but has been introduced to reduce the 
grid computations. Our best-fitting model has a log-likelihood 
£„, = -314.17. The 95% confidence intervals are determined 
by all models whose log-likelihood £ is 2(£„ - £) < ^'15(0'), i-e. 
2({,„ -{) < 24.99579 or -326.6 < £ < -314.17. If we want 
to test if a model with 3 episodes of high mass loss yields a 
significantly better result than with 2 high-mass epochs, we are 
introducing 4 extra pairs (/?,, M,), i.e. 8 free parameters. Hence, 
the more elaborate model with 23 parameters is preferred when 
the maximum of the log-likelihood of the 3-shell model {£f^) 
obeys 2(C - O > xl(») or f^f > -306.41. 

4. Results 

4.1. Modeling the CO rotational line profiles 

Comparing the observed CO line profiles with theoretical pro- 
files predicted by a model assumin g a constant mass loss con- 
firms the hypothesis proposed by iKemper et al.l (l2003h that 
variable mass loss may play an important role in the shaping 
of the spectral line profiles (see dotted line in Fig.|3]l. 

Using the procedure as outlined in Sect 13.1.31 a total 
amount of ~300 000 models was run to estimate the mass-loss 
history of WX Psc. A luminosity of 1 x IO'^Lq was assumed, 
which together with the integrated luminosity of the SED yields 
a distance of 833 pc. The expansion velocity traced by all line 
profiles is = 18 + Ikms"'. For the carbon and oxygen 
abundances in the outer atmosp here of the target, we opt to use 
the cosmic abundance values of lAnders & Grevessd (1 19891 see 
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Table [3]l, and it is assumed that in case of the oxygen-rich tar- 
get WX Psc all of the available carbon will be locked in ccB 
Other model parameters as stellar temperature (or stellar ra- 
dius Ri,), dust condensation radius /?inner, outer radius of the 
envelope Router and the mass-loss profile were kept free in the 
grid running. In case of constant mass loss, the value of /Pouter 
is set at the radius where the CO abundance drops to 1 % of 
its value at the photosphere, using the photodissociation results 
oflMamonet al. (1988) (see Decin et al. 2006 ). How ever, the 
CO photodissociation radius by iMamon et alJ (il988h were ob- 
tained in case of constant mass loss, and it is questionable how 
this applies when dealing with variable mass loss. In case of 
variations in the mass loss, we therefore opt to have Romer as 
free parameter, with the additional constraint that Romei should 
be smaller than the /?outer-value obtained using the formula of 



' Note that in case of a varia ble mass loss, t he [C O/H2]-ratio is as- 
sumed to follow the results of iMamon et al.l jl988l) . Since the total 
abundance of the main molecule H2 varies according to the M-profile, 
the total CO abundance varies accordingly. 



iMamon et al. ( 1988 ) for a value of the mass loss equal to the 
mass-loss rate at the outermost radius point. We hereby want to 
ensure not to go beyong the photodissociation radius as calcu- 
lated bv iMamorTetaLl IT988I) . Moreover, it is always checked 
that the full lin e forming region ( as traced by I{p) x p^, see e.g. 



Fig. 5 and 6 in lDecin et al. Il200d) of all rotational CO lines lies 

within /Pouter- 

Parameters characterizing the CSE of the model with the 
best goodness-of-fit are specified in the second column of 
Table[3] The derived temperature profile, velocity structure and 
mass-loss history are displayed in Fig. |2j a comparison be- 
tween observed and theoretical line profiles is shown in Fig. [3] 
We note that the abruptness of changes in the mass-loss profile 
causes the small-scale ripples seen in the CO line profileo The 
presence of the noise in the line profile does however not influ- 
ence the conclusions in this paper. The dust-to-gas ratio is de- 



- Theoretical profiles are calculated with 150 impact paramters and 
200 quadrature point in the /y/? dp-integral. 
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Table 3. Parameters of the model with best goodness-of-fit for 
WX Psc. In the second column models parameters as derived 
from the CO rotational lines are listed, in the third column the 
parameters as estimated from the SED fitting are given. The 
numbers in italics are input parameters that have been kept 
fixed at the given values. There is, of course, an interdepen- 
dence between T*, and L* and between L*, the distance, 
and nihoi from which the distance has been derived. 



parameter 


from CO lines 


from SED fitting 


[K] 


2000 


2000 


R^, [10" cm] 


5.5 


5.5 


L. [10*Lo] 


1 


1 


dust-to-gas ratio 


0.004 


0.004 


£(C) [10-*] 


3 




£(0) [10-'*] 


8.5 




distance [pc] 


833 


833 


dinner [^*] 


5 


5 


^outer [R*] 


1750 


1500+ 


Voo [kms"'] 


18 


18 


M(r) [Moyr-'l 


(A)~610-'5 (-5-30 i?J 


~4 10-5 (^5_50i?j 



(6-5) line is somewhat too high. This may indicate that the 
temperature in region C, where the dominant part of these lines 
is formed, may be somewhat lower than is inferred. 

Almost half of the models was run with a third shell ex- 
tending between ~100 and ~600/?*, and with mass-loss rates 
ranging between 2 x 10"^ and 3 x 10"^ Mq yr"' . The model with 
the best fit had M -5 x 10"^ Mq yr"' both in region (B) and in 
this extra third shell, i.e. the statistical analysis gives preference 
to a 2-shell model above a 3-shell model. Models still having an 
acceptable fit to the line profiles, all had M < 2.5x10-^ Mq yr"' 
in this third shell. 



(B)~510-*' (-30-850 



4.2. SED fitting 

The SED of the best fitting model is presented in Fig. [T] and 
the implied model parameters are given in Table [3] As de- 
scribed above, the model consists of two shells that present two 
episodes of high mass loss. The intermediate period of signifi- 

(C) ~5 10-5 (-850-1450 i?*) -3 10-" (-800- 1500 cantly lower mass loss is not modeled. In general, we find that 

(D) -3 10 (-1450-1750/?,)- j-jjg y/ijid structure that best reproduces the observed SED has 

very similar parameters to those derived from modeling the CO 
rotational lines. In particular, the two methods agree on the lo- 
cation of the inner edges of the two shells. The outer radii of 
the shells can be chosen to be roughly consistent with the val- 
ues derived from the molecular modeling. It is impossible to 
derive statistical uncertainties on the derived parameters since 
the best fitted model is not obtained through a formal fitting 
procedure (like was done for the CO lines). Such a formal fit- 
ting procedure is inhibited by the very unequal weighting that 
is given to the various regions of the SED. For example, the re- 
gion which is dominated by water vapour emission and absorp- 
tion is virtually neglected. While the integrated photospheric 
component is used as a constraint, no single optical wavelength 
range is given priority. In contrast the 8-25 fim range, which 
exhibits the clearest dust spectral signatures, is inspected in de- 
tail. Furthermore, the IR-submm slope is used as a strong con- 
straint eventhough it is determined by far fewer points than the 
mid-IR part. However we can determine by how much we can 
vary the parameters without significantly reducing the quality 
if the fit. Varying, e.g., the inner radii of the two shell model 
by ~10 per cent, significantly degrades the quality of the fit. 
We also note that the dust model results are less sensitive to 
the precise extent of the dust shell than to the amount of mass 
contained in the shell. Therefore, the outer radii are less well 
constrained. An important conclusion that can be drawn is that 
the self-absorption and relative strengths of the 10 and 20 /zm 
silicate features arise from a dense and compact shell close to 
the photosphere. To illustrate this point we also show in Fig.[T] 
the "best" fitting model using a single shell from a constant 
mass loss (M- 1.5 10-^ Mq; 7;=5-200 Rq). A more extended 
shell would yield too much cold dust and a less dense wind 
would not yield enough self-absorption. This compact shell is 
however not sufficient to explain the far-IR data. The slope of 
the SED in the far-IR is dominated by cool material (see inset 
in Fig.[T]i. 



' The outer radius mentioned in the SED fitting refers to the outer 
edge of the second shell, called region (C) in Fig. |2l and should be 
compared to the 1450 i?* value derived from the CO modeling. 



rived from the mass-loss rate and terminal velocity and equals 
0.004. This val ue for the dust-to- gas ratio is quite uncertain. 
As described in Decin et al. I (l2006l) . this value is taken constant 
throughout the entire envelope. However, a higher mass-loss 
rate as in region (C) can drive a dusty wind with a different 
velocity than in a low-mass region and perhaps also with a dif- 
ferent dust-to-gas ratio. In this respect, we note that the differ- 
ent CO lines, originating from various parts in the wind, do not 
exhibit different line widths (Fig.[3]l. 

Before comparing with the results of the SED fitting 
(Sect. 14.21 ). we first discuss the derived M-profile and associ- 
ated uncertainties (see bottom panel in Fig.|2]l. Acceptable val- 
ues for the mass-loss rate in region (A) range from 5 x 10-^ 
to 1 X lO-^Moyr-', for region (B) from 2.5 x 10-^ to 8 x 
10-** Mo yr-\ for region (C) from4x 10-^ to 6.5x 10-^ Mg yr-', 
and for region (D) from 1 x 1 0-^ to 4 x 1 0-* Mq yr- ' . Especially 
region (C) with enhanced mass loss is constrained quite well 
since the (high-quality) low-excitation lines are mainly formed 
in this region. Notice that the temperature increase in region 
D (see top panel in Fig. |2]| is quite uncertain: an altered tem- 
perature structure with a continuously decreasing temperature 
reaching lOK at the outer radius, yields line intensity profiles 
which are better for all low-excitation transitions, and has a log- 
Ukelihood that is significantly better than that of our 'best-fit' 
model with an amount Af = 105. The main cause of uncer- 
tainty in T(r) is that interactions between the different shells 
and with the ISM are not taken into account. 

Inspecting Fig. [3] we see a good resemblance between the 
observations and model predictions for the CO (1-0), (2-1) and 
(3-2) Unes. The predicted peak intensity of the CO (4-3) and 
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4.3. CO versus SED fitting 



5. Discussion 



The fact that these very different tracers yield such similar pa- 
rameters with respect to the location of the densest shells is 
remarkable and builds further confidence in the reality of the 
mass-loss variabilities. However, it should be noted that the 
mass-loss rates derived from fitting the SED (column 3 in 
Table [3]) are higher than those found from modeling the CO 
emission (column 2 in Table |3]l, by roughly a factor 6. Partly, 
the discrepancy can be accounted for by considering the ma- 
terial contained in the regions that are not included in the dust 
model (B and D). Some experimenting with the dust model- 
ing shows that the differences can not be reconciled by adapt- 
ing the dust composition alone. Upon perusal of Table [1] it 
appears that those studies that take the thermal IR emission 
into account systematically derive a denser wind than those 
using gas tracers alone. The difference may be related to the 
assumed CO/H2 ratio (see TableO. Under the assumption that 
all carbon is locked in CO, the CO/H 2-ratio used in this paper 



is 3 lO"'*. IZuckerman & Dvckl (Il986l) used a value of 5 10 ^^ 



for the CO/H2-ratio in O-rich giants, while iKnapp & Morris 
(119851) reported a value of 3 10"^. This value differs from AGB 
star to AGB star, as it is dependent on the amount of dredge- 
up. We estimate the uncertainty on the used CO/H2-ratio to be a 
factor 2. In addition, as described in Sect. 14. II the derived dust- 
to-gas ratio as calculated from the CO modeling is somewhat 
uncertain. This value is not only taken constant throughout the 
whole envelope, but moreover it is assumed that all the dust 
particles participate in accelerating the wind to the observed 
outflow velocity. The discrepancy in dust content between the 
two types of models may signify that a sizable fraction of the 
dust, though present, does not contribute much to the accelera- 
tion of the wind. If some fraction of the dust particles is effec- 
tively shielded from illumination by the star, they will not feel 
the radiation force from the star and, by consequence, will not 
contribute to the acceleration of the wind. This kind of shield- 
ing can typically occur in randomly distributed density inho- 
mogeneities (usually referred to as clumps) or can be due to 
deviations of spherical symmetry. Only some fraction of the 
dust is hence responsible for the driving; via collisions with the 
inert gas, the gas is accelerated, which in it turn may acceler- 
ate the inert dust. In the acceleration of the gas particles, it is 
not only assumed that all dust particles participate in this pro- 
ces, but also that the dust grains transfer all of the momentum 
they acquire from the radiation field to the gas through colli- 
sions. Therefore, an alternative explanation could be that there 
is an incomplete mome ntum coupling between the dust and 
the gas in the wind (e.g. MacGregor & Stencell 1992 ). Lastly, 
an explanation could have been dust which is not partaking in 
the outflow but present in a stable configuration, i.e. a disk. 
However, the gaseous material present in such a slowly rotat- 
ing disk would be detected via the presence of a very narrow 
line profile (typically line width less than 5kms ') superim- 
pose d on a broad outflow component (e.g. as is the case of RV 
Boo, Bergman et al. 2000l) . Since such a narrow component is 
not detected in any of the line profiles this explanation is ruled 
out. 



The diagnostic power of estimating the mass-loss history from 
modeling the rotational CO Une profiles and the full SED, is 
evident. Both our CO and SED modeling are however based 
on one common assumption, namely a spherically symmet- 
ric envelope. We comment on this assumption in Sect. 15.11 In 
Sect. 15.21 we compare our results to other studies. We discuss 
the time scale of the mass-loss variability and its implication 
on the possible mechanism driving this variability in Sect. 15. 31 

5.1. Structure of the circumstellar envelope 

The detailed structure of the CSEs of AGB stars and the physics 
driving this structure remain unknown after half a century 
of study. The main reason being that the cool material is a 
real challenge to high-resolution observations. Fror n a survey 



of AG B envelopes using millimeter interferometry, iNeri et al 



(119981) concluded that most (~70 %) AGB envelopes are con- 
sistent with spherical symmetry. WX Psc belongs to the small 
group of AGB stars whose CSE has already been studied us- 
ing different techniques, and for which several indications are 
fou nd for dev i ations from a spherically symmetric envelope. 
(^7) lNeriet alld 1998b fitted the CO (1-0) and (2-1) interfero- 
metric observations with a two-component Gaussian visibility 
profile consisting of a circular and an elliptical component. The 
circular flux-dominant component corresponds to a spherical 
envelope with a 29". 6 diameter, the elliptical one to an ax- 
isymmetric envelope with major axis of 9". 8 and minor axis 
of 6". 8, and a position angle of -45° (E to N). (2) Broad 



band n ear-infrared polarimetry performed bv lMcCall & Hough 



(Il980l) revealed that WX Psc is highly polarized in the I and 
/-band, while the polarization pattern is not dominant in the 
H and ^T'-band. This may either indicate an alignment of dust 
particles or an asymmetrical dust-shell structure. (3) Another 
technique of imaging the circumstellar dust is f rom bispectrum 



speckl e-interferometry. Using this technique, iHofmann et al 
( 200ll) showed that the H and /T'-band images of WX Psc ap- 



pear almost spherically symmetric, while the 7-band shows a 
clear asymmetry. Two structures can be identified in the /-band 
image: a compact elliptical core and a fainter fan-like structure, 
along a symmetry axis of position ang l e -28° , out to distances 
of ~200mas. (4) iMauron & Huggins (l2006h imaged the cir- 
cumstellar dust in scattered light at optical wavelengths. While 
the images of the extended envelope appear approximately cir- 
cularly symmetric in the V-band, the Hubble Space Telescope 
(HST) images revealed that the inner region of the CSE consists 
of a faint circular arc at ~0.8"to the NW and a highly asym- 
metric core structure, with a bright extension out to ~0.4"(or 
90 at 833 pc) at position angle -45°. Most likely, the HST 
image captures the extension of the core asymmetry seen in the 
7-band. They concluded that WX Psc is somewhat similar to 
IRC -1-10216 where approximate circular symmetry and shells 
in the extended CSE co-exist with a strong axial symmetry close 
to the star. 

Nevertheless, we assume a spherically symmetric geometry 
for the whole wind structure when modeling the SED and CO 
line profiles. The different direct pieces of evidence for devia- 



L. Decin et al.: The variable mass-loss of the AGB star WX Psc 



9 



tion of spherical symmetry all trace small scales 200 mas). 
On larger scales, the shell does not appear to deviate signif- 
icantly from spherical symmetry. Note that the beam profiles 
used in the CO observations are larger than 6", and thus mea- 
sure the temperature, density and velocity averaged in all direc- 
tions in the regions where the lines are mainly formed. Hence, 
the main stellar and mass-loss properties can be inferred from 
spherically symmetric models in fair approximation. In reality, 
these average properties are formed from a superposition of all 
these structures in this region. 

5.2. Comparison with other studies 



In Sect. 15.2.11 we compare our results with mass-loss parame- 
ters obtained in previous studies also tracing the mass-loss vari- 
abihty of WX Psc. In Sect.l5.2.2l we focu s on t he V, J, H, and 



K' im ages of iMauron & HugginsI (12006b and iHofmann et al 



( 120011) probing both the dust-scattered stellar and galactic Ught 
and the thermal dust emission. 



5.2.1. Mass-loss parameters 

Hofmann et al. I (1200 ih modeled the SED between 1 jum and 



1 mm, together with the visibility function s at 1.24 (7 band) , 
1.65 (H), 2.12 (K') and llyifm using dusty dlvezic et al.ll 19971) . 
The visibility data provide information on the innermost re- 
gions of the circumstellar environment, probing the scattering 
properties of the dust (at J) and thermal emission properties 
of the hot and warm dust (at H, K and 1 1 //m). Assuming a 
spherical outflow, they found that both the near-IR visibilities 
as well as the mid-IR photometric points up to lOOyum could 
be well reproduced by a two-component outflow in which 
dust condenses at 900 K. The most recent mass loss, charac- 
terized by a constant velocity outflow (i.e. p oc r"-), is M = 
1.3 - 2.1 X 10"'' Moyr"'. At distances beyond 135/?* the den- 
sity decrease is shallower, matching p ~ r (see the dashed 
black line in Fig.|4]i. At wavelengths beyond about 400 yum this 
model underestimates the observed flux. 

The y-band data reveal an asymmetry, and the 7/-band vis- 
ibility shows a puzzling drop in the //-band at spatial fre- 
quencies ^ 14 cycles per arcsec, corresponding to structure 
smaller than the cond ensation radius. The spherical model of 
Hofmann et alJ (1200 ih can not address these observa tions. The 
above issues were taken up by lVinkovic et al.l (120041) . They ex- 
plain the image asymmetries as originating from about 10 Mq 
of wind material that is swept-up in an elongated cocoon whose 
expansion is driven by bipolar jets and which extends out to a 
radial distance of ~ 1 100 AU (or 300 /?*), with an opening angle 
of ~30°. The two cocoons have a characteristic lifetime that is 
^ 200 yr. To explain the mid-IR spectrum they invoke a spher- 
ical envelope (in which the cocoons are embedded) containing 
-0.13 Mq of material extending out to -35" (or -8 000/?*). 
The size of the envelope corresponds to a dynamical flow 
time of 5 500 yr. Over time the mass loss rate has decreased. 
While M ~ 2 X 10"^ Mgyr"' until 550 yr ago, the most re- 
cent mass-loss is estimated to be only -9x10"'' yr"' . Two- 
dimensional radiative transfer calculations demonstrate that 



this model can explain both the 7-band asymmetry, being dom- 
inated by scattered light escaping through the cone, and the 
spherical symmetry seen in the H and K'-bwd images, being 
due to predominantly thermal dust emission. This model also 
fits the general structure of the SED quite well, but has some 
problems in explaining the silicate feature around 10 yum and 
the photometric data at /i > 100 ;um. 

In ge neral, the mass- l oss rat e in t he spherical envelope de- 
rived by ' Hofmann et"al] (I2OOI') and 'Vinko vic etall (|2004 is 
much larger than what is implied by the CO lines (see Fig. |4|l 
and is not consistent with the observed CO rotational line pro- 
files (see Fig. Is). The CO (1-0) through (7-6) lines are a sen- 
sitive probe of a geometrically extended part of the CSE, save 
for the innermost region. This again may indicate that part of 
the dust does not participate in driving the wind. 
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Figure 4. Number densities as a function of radial distance 
from the star as deduced by different authors. Full black 
line: spherical CSE of this study; dashed black line: spher- 
ical CS E of Hofmann et al. (2001); full grey line: spherical 
C SE of IVinkovic et all (|2004 : dashed grey line: bipolar cone 
of Vinkovic et a l.' ('20041: dotted grey line: number density of 
bipolar cone of iVinkovic et al., (.2004) scaled to a sphere. 



5.2.2. Dust scattered light 



Beside s the work by iHofmann et al. ( 2001 ) and Vinkovic et al 



(|2004 . HST observations bv iMauron & HugginsI ( I2OO6I) also 
trace mass-loss variability in the ambient environment of WX 
Psc. The HST-image shows part of a faint circular arc at 
-0.8"(or -ISO/?*) from the center. We performed additional 
simulations introducing a third shell, positioned at this loca- 
tion (see Sect. 14.1b . These simulations demonstrate that it may 
be possible that a faint arc is present with a mass-loss rate 
^ 2.5 X 10"^ Mq yr"', although the log-likelihood of this more 
sophisticated model is not significantly better than for the sim- 
pler two-shell model. The presently available rotational CO 
line profiles can hence not be used to draw any conclusions 
on the presence of this faint arc. 
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Figure 5. CO rotation al line profiles of 
WX Psc (full grey l ine; iNvman et al 111 9921: 
lKemperetai]|2003l) compared with model 
predictions based on the param eters as 
proposed by Hofmann et al. ( 2001) (dotted 
black line) and .Vinkovic et aL (.2004) (full 
black line). 



Mauron & Hugginsl(l200 6') have obtained the B and V-band 
scattered light profiles of WX Psc. The scattered light is pre- 
dominantly due to external illumination of the shell. These im- 
ages show a relatively smooth decrease of intensity as a func- 
tion of distance and at face-value do not seem to necessitate the 
invocation of mass-loss variability to explain them. However, 
the observed scattered light profile is very sensitive to the loca- 
tion of the external illumination sources and the phase function 
of scatterers. Formally the expression that these authors have 
used to fit to the observed profiles under the assumption of uni- 
form illuminations is incorrect. The expression that they ap- 
ply concerns the scattered part of the observed intensity, while 
under uniform illumination this part exactly compensates the 
intensity which is scattered out of the beam. It is at present 
unclear how a variable mass-loss history would translate into 
an observed intensity profile on the sky as a function of these 
parameters. It is obvious that a forward peaked scattering func- 
tion will tend to de-emphasise any mass-loss variations at large 
distances. But whether this effect is sufficient to explain the 
smooth brightness distribution remains to be seen. A proper 
model for the light scattering in the CSE of WX Psc, taking 
into account its position in the galaxy, should be constructed. 
This, however, is clearly beyond the scope of this paper. 

5.2.3. Conclusion 

To conclude, we use constraints from the dust and the CO 
molecules in the shell to determine the density profile within 
the few to few thousand stellar radii regime. The tracers we use 
are sensitive to these scales in a series of overlapping distance 
steps. The real power of our modeling comes from using both 
the spectral information and the absolute strength of the ob- 
served quantities, be-it the molecular lines or the dust excess. 
These are all explained to satisfactory detail by the model pre- 



sented here. The limitations of the method are also clear. While 
this yields good estimates for the average density as a function 
of distance, neither the large-beam observations nor the spheri- 
cal model have the power to draw conclusion about small-scale 
{AR/R X 0.5) structures or asymmetries. 

5.3. Mass-loss variability 



As already noted in Sect. 15.11 auron & HugginsI (l2006h con- 
cluded that WX Psc resembles IRC H- 10216 where approximate 
circular symmetry and discrete shells exist in the extended en- 
velope. In case of IRC -1-10216, the shell spacing varies be- 
tween 200 and 800 yr, a lthough intervals as short as 40 yr are 
seen close to the star (^ Mauron & Hug gins 2000l) . The time 
spacing between shell (A) and (C) in WX Psc is ~800yr. The 
total envelope mass is estimated to be 0.03 Mq, of which ~95 % 
is located in region (C). Assuming WX Psc to be a typical 
(solar-type) AGB star losing some 0.5 Mq during its AGB evo- 
lution lasting some 100 000 yr, WX Psc has already lost ~6 % 
of the total mass that it will shed in its AGB-phase in a period 
of 600 yr (region C). Such a strong ejection event may indi- 
cate that WX Psc is currently in its superwind phase. Assuming 
that a typical oxygen-rich AGB star expels 0.1 Mq during the 
last lOOOyr that it expe riences a superwind mass-loss phase 
dMarigo & Girardill2007h . few such high mass-loss phases may 
occur during the superwind phase of WX Psc. 

Various models have been proposed for the origin of dis- 
crete shells in CSEs. One type of mod els is based on the 



effects of a binary companion (see, e.g., iHarpaz et al. 1997 



Mastrodemos & Morris..l999). A different scenario is invoked 



by ISimis et al.l (120011) . who suggested that such shells are 
formed by a hydrodynamical oscillation due to instabilities 
in the gas-dust cou pling in the CSE while the star is on the 



AGB. ISoked ( l2000l) invoked another mechanism, being a solar- 
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like m agnetic activity cycle in the progenitor AGB star. ISoker 
(l2006h raised the possibility that semi-periodic oscillations in 
the photospheric molecular opacity may be another candidate 
mechanism for the formation of multiple semi-periodic arcs. 



6. Conclusions 

In this paper, we have analyzed the rotational CO (1-0) 
through (7-6) line profiles of WX Psc in combination with the 
dust characteristics in the full SED from 0.7 to 1300jum. This 
combined analysis yields strong constraints on the density and 
temperature structure over a geometrically large extent of the 
circumstellar envelope, save for the innermost region. Both CO 
line fitting and dust modeling evidence strong mass-loss modu- 
lations during the last 1700 yr Particularly, WX Psc underwent 
a high mass-loss phase {M~ 5 10"^ Mo/yr) lasting some 600 yr 
and ending ~800 yr ago. This period of high mass loss was fol- 
lowed by a long period of low mass loss (M~ 5 10"** Mo/yr). 
The current mass loss is estimated to be ~ 6 10"^ Mo/yr. We 
want to emphasize that the time spacing of these mass loss 
modulations are independently well constrained by both the 
molecular gas and dust modeling. The uncertainty on the mass- 
loss rate in the different phases is estimated to be a factor of a 
few, mainly caused by the uncertainty on the dust-to-gas ra- 
tio and the CO/H2 ratio, and by the assumption of a homoge- 
neous spherically symmetric envelope. Most importantly, the 
mass-loss history that we derive is significantly different from 
a constant mass-loss rate model, and moreover implies that WX 
Psc can have at maximum in the order of a few high mass-loss 
phases during its final superwind evolution on the AGB. 
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